Coulomb Systems with Ideal Dielectric Boundaries: 
Free Fermion Point and Universality 

B. Jancovici 1 and L. Samaj 1 ' 2 
February 1, 2008 

Abstract 

A two-component Coulomb gas confined by walls made of ideal dielectric ma- 
terial is considered. In two dimensions at the special inverse temperature (3 = 2, 
by using the Pfaffian method, the system is mapped onto a four-component Fermi 
field theory with specific boundary conditions. The exact solution is presented for 
a semi-infinite geometry of the dielectric wall (the density profiles, the correlation 
functions) and for the strip geometry (the surface tension, a finite-size correction of 
the grand potential). The universal finite-size correction of the grand potential is 
shown to be a consequence of the good screening properties, and its generalization 
is derived for the conducting Coulomb gas confined in a slab of arbitrary dimension 
> 2 at any temperature. 
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1 Introduction 



Very recently, the interface between a Coulomb gas and a dielectric wall with zero dielectric 
constant (termed ideal dielectric) has been studied by one of us (L.S.) |l| in the case of an 
exactly solvable two-dimensional model 0: a two-component plasma of point-particles. 
The surface tension has been obtained, for an arbitrary temperature within the range of 
surface stability of the model (3 < 3 ((3 is the inverse temperature). This range exceeds 
the bulk range of stability j3 < 2. However, the method used in a mapping onto a 
solvable boundary sine- Gordon field theory, provided only the surface tension but did not 
give any information about the microscopic structure of the interface: density profiles and 
correlation functions. 

It can be expected that this microscopic information is obtainable at the special in- 
verse temperature (3 = 2. This is because, at this special temperature, for solving other 
problems, the model can be mapped onto a free fermion field theory (This temperature is 
the boundary of the stability domain of the model: for point particles, at a given fugac- 
ity, the bulk density and other bulk thermodynamic quantities diverge, however they can 
be made finite by the introduction of a small hard core in the interaction. Anyhow the 
truncated many-body densities are finite, even for point-particles). The two-dimensional 
two-component plasma dX (3 = 2 has been extensively studied, in the bulk ||, and also 
near a variety of possible interfaces, in particular the interface with a plain (simple) hard 
wall |4j] and with an ideal conductor wall - • However, the case of the interface with 
an ideal dielectric could not be solved by the simplest possible extension of the mapping 
on free fermions with an inhomogeneous "mass" term which had been used for other 
interfaces. 

In the present paper, we do solve the problem of the interface with an ideal dielectric 
by an appropriate generalization of the mapping. In the previously studied cases, the 
mapping was obtained by associating a two-component Fermi field with each point of 
space. In the present case, inspired by a remotely related reference ||, we realized that a 
mapping on a four-component Fermi field is needed. Through that mapping, we are able 
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to compute the density profile and the correlation functions near the interface. This is 
described in Section 2 dealing with a semi-infinite geometry of the dielectric wall. 

In Section 3, we turn to the problem of the two-dimensional model at /3 = 2 confined 
in a strip of width W between two ideal dielectric walls. We show that the grand potential 
has the same universal finite-size correction of order 1/W as in the previously studied case 
of ideal conductor walls @. Actually, as shown in Section 4, this finite-size correction is 
a consequence of the good screening properties and its generalization can be derived for 
a conducting Coulomb system of arbitrary dimension at any temperature. 

2 Model and method of solution 
2.1 Pfaffian method 

We consider an infinite 2D space of points r e R 2 defined by Cartesian (x, y) or complex 
[z — x + iy, z — x — iy) coordinates. The model dielectric-electrolyte interface is localized 
along the x axis at y — 0. The half-space y < is assumed to be occupied by an ideal 
dielectric wall of dielectric constant e = 0, impenetrable to particles. The electrolyte 
in the complementary half-space y > is modeled by the classical 2D two-component 
plasma of point particles {j} of charge {qj = ±1}, immersed in a homogeneous medium 
of dielectric constant = 1 and interacting via Coulomb interaction. 

At the special inverse temperature /3 = 2, in order to avoid the collapse of positive 
and negative particles, we start with a lattice version of the Coulomb model. There are 
two interwoven sublattices of sites u & U and v G V. The positive (negative) particles sit 
on the sublattice U (V); each site can be either empty or occupied by one particle. We 
work in the grand canonical ensemble, with position-dependent fugacities ((u) and £(i>) 
in order to generate multi-particle densities. In infinite space, the Coulomb potential Vq, 
at spatial position r, induced by a charge at the origin, is vo(r) = — ln(|r|/a), where the 
length constant a fixes the zero point of energy. The presence of a dielectric wall induces 
to a particle at r = (x, y), or z, an image at r* = (x, —y), or z, of the same charge 0. We 
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consider only neutral configurations of iV positive and N negative particles, the complex 
coordinates of which are itj and v^, respectively. The corresponding Boltzmann interaction 
factor, involving direct particle-particle interactions qiqjV$(\ri — rj\) and interactions of 
particles with the images of other particles qiC[jVo(\ri — as well as their self-images 
(l/2)g^ (|n-r*|),is, at /3 = 2, 



B N ({u t }; { Vi }) = (-lfa 2N H( Ui - UiXvi - v t ) 

(2.1) 



Ui<j(Ui ~ Uj)(Ui - Uj)(Ui - Uj)(Ui - Uj)(Vi - Vj)(Vi - Vj)(Vi - Vj){Vi -Vj) t 



Ui,j(Ui - Vj)(Ui - Vj)(Ui - Vj)(Ui - Vj) 

By using an identity of Cauchy JTIJ 

det ( 1 \ = ( _ 1)W n^-f)(j-4) (2.2a) 

V' 4 * 3 / i,j=l,...,2N Uij'Vi 

with the identification 

^2i-i = Mj, z 2 j = Mi, = u i, 4i = ^i, « = 1, • • • , N (2.2b) 



the Boltzmann factor (|2.1| ) is expressible as the determinant of an N x N matrix whose 
elements are 2x2 matrices: 



B N ({ui}; {vi}) = det 



The grand partition function reads 



/ a a \ 

Ui — Vj Ui — Vj 



(2.3) 

V U % - Vj Ui - Vj J N 



E = l+J2C(u)C(v)B 1 (u;v)+ J2 ((ui)((u 2 )((v 1 )((v 2 )B 2 (u h u 2 ,v h v 2 ) + ... (2.4) 



where the sums are defined in such a way that configurations which differ by a permutation 
of identical-charge particles are counted only once. 

We introduce a couple of Grassmann variables for every site u G U, (if)^, ■0^), and every 
site v G V, (ipl, if)%), and order all of them into a single vector t Q = (. . . i/)^, ij)^ . . . . . .) 

whose components obey the ordinary anticommutation rules {6^, 6j} = 0. Let us consider 
the Grassmann integral 

Z = J d9 e~ So (2.5a) 



where dO = U v d^d^ U u dipu^u and 



-S = - "OAO 
2 



(2.5b) 



is (minus) the Gaussian action with an antisymmetric matrix A, Aij + Aji = 0. We recall 
that A -1 is also antisymmetric. The averaging with So will be denoted by 
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(2.6) 



Gaussian integrals of type ( |2.5|) are expressible as JTT 

Z = Pf (A) 



(2.7) 



where Pf(A) is the Pfaffian of the antisymmetric matrix A, satisfying the well known 
identity 

Pf(A) 2 = detA (2.8) 
The two-variable averages are given by 



(2.9) 



The standard Wick theorem for fermions generalizes to higher-order averages of 9 variables 



12fl . The A-matrix will be chosen such that its inverse satisfies the following relations: 

(A 



-l\aa' 
Juu' 



(CC) = o 



(2.10a) 
(2.10b) 



i.e., the sites of the same sublattice-type do not interact with each other, and 
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(2.10c) 



(2.10d) 



with the required antisymmetry property (A : )"^ = — (A Here, the rows and the 

columns of the matrices are numbered as a, (3 = 1, 2. 



We now introduce an antisymmetric "mass" matrix, 



M aae 



M af3 
uv 

MS? 








1 o K(u) ^ 



-K{u) o 



/ 



' o i((v) ^ 

v -iC(«) J 



(2.11a) 

(2.11b) 

(2.11c) 
(2.11d) 



Define 



with the action 



/ 



dd e" 



-S = i*#(A + M)# 



l./,2 
v 



(2.12a) 



(2.12b) 



Clearly, 

Z = Pf(A + M) (2.13) 

Let us expand the ratio Z/Z Q in fugacities {((u)} and {((v)}. Owing to ( |2.10| ), only 
neutral contributions with an equal number of {C(w)} an d {C( v )} wm survive, 



A? 



A' 



y co 

y = i + e (-I) 7 " e n c(«oc(«o ( n («« 

^0 N=\ U!<u 2 <— <u N €U i=\ \i=l 



;2.i4) 



With the definition of pair averages ( |2.10|) , the Wick theorem for fermions (see, e.g, fll 
T"2|j ) implies 



/ N 

n {«« 



N 



\i=l 



\i=l 



n {«« 



( 



-l^det 



;2.i5) 



V Ui - Vj u % - Vj J iJ=1 
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Comparing with ( |2.3|) and (|2.4| ) one concludes that 

„ Z Pf (A + M) 



Pf(A) 



(2.16) 



Introducing the matrix K = MA 1 , ( |2.16D can be rewritten as 



In 5= -lndet(l + K) = ^Trln(l + K) 



(2.17) 



Marking the charge sign of the particle at r = (z,z) by an index s = ±, i.e., identifying 
u = (r, +), v = (r, — ) and £(«) = C+( r )? C( v ) — C-( r ); the K-matrix has the elements 



/ a 



a \ 

Ktf(r, r') = iC.(r) ^/ ^ 

V 2; — 2' z — z' / 
where the a, a' = 1, 2 indices label the elements of the 2x2 matrix. 



(2.18) 



2.2 Many-particle densities 



The calculation of many-body densities from the generator of type ( [2.17D is straightforward 
(see, e.g., When S is the area of a lattice cell, using the equality 





(2.19) 



one finds for the density of particles of one sign (number of such particles per unit area) 



1 d 

5 «Ci(n) 
z Oil 



Ur)=C 



Here, the matrix 



G 



K 



(2.20) 



(2.21) 



27ia( 1 + K 

is defined with ^-elements evaluated at constant fugacity, ( s (r) = (, and the rescaled 
fugacity m is defined by m = 2na(/S. Using 

d 1 1 dK 1 



9G 2 (r 2 )l+K 1 + K 9G 2 (r 2 ) 1 + K 
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[2.22) 



one finds for the truncated two-body density 

1 



If 



in: 



Ur)=( 



cti,a 2 

By successive iterations, one finds for the truncated fc-body density 



(2.23) 



< )T Jn,..,^ = H) fc+1 ™ fc E E G^{ ni ,r i2 )...G^{r ih ,r k 

a u ...,a k ( n i 2 ...i fc ) 



(2.24) 



where the second summation runs over all cycles {i\i 2 ■ ■ - ik) built with {1, 2, . . . , k}. 

In the continuum limit where the lattice spacing goes to zero, the G ss > (s, s' = ±) 
matrices, defined by (|2.21|) , satisfy the integral equations 

I -L l —\ 



Gssin, r 2 ) 



G- a>8 (r\, r 2 ) = i 



i (-) I dV 



2tt/ JD 

{ 1 



Z\ — Z Z\ — Z 
-1 -1 



V 



G_ s , s (r,r 2 ) (2.25a) 



1 

2^ 



Z\ — Z Z\ — Z ' 

-J—\ 



Z\ — Z 2 Z\ — z 2 

-1 -1 



-i (-) I d 2 r 



\ Z\ — Z 2 Z\ — z 2 

(-L- 



V27T/ 



D 



\ 



Z\ — Z Z\ — z 

-1 -1 



G ss (r,r 2 ) (2.25b) 



V Z\ — Z Z\ — z I 

where the domain of integration D is the half plane y > 0. The operators 



d 



1 d_ _ .d_ 

2 I dx dy 



d. 



1 /_£?_ . jT 

2 I <9x dy 



act as follows [11 



d z 



1 



- = 7T<5(r — r'), 



i 



z — z z — z 

Differentiating appropriately equations (|2.25| ), one gets 



- = 7i5(r — r') 



( 



d 



21 



-d Zl 



o a 



21 



-3* 



My) G - S)S ( r i, r 2 ) 



1 



G_ a> .(n,^) = -i-5(n-r 2 )l + i - G ss (n,r 2 ) (2.28b) 



/ 



??? 



(2.26) 



(2.27) 



(2.28a) 
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where 1 denotes the 2x2 unit matrix. Using d z d s = (1/4) A, these equations can be 
combined into 

(_Aj + m 2 )G ss (r 1; r 2 ) = rn5{n - r 2 )l (2.29) 
The boundary conditions are given by the integral equations ( 12.251) : 



C(n,r 2 ) + C(n,r 2 ) = 0, 

Gl 2 sl ( n ,r 2 ) + GfAn,r 2 ) = 0, 



[n V r 2 ) G dD 
[n V r 2 ) G &D 



(2.30a) 
(2.30b) 



where dD is the domain boundary y = 0. 

Since G ss '(ri, r 2 ) is translationally invariant along the x axis, for solving the partial 
differential equations ( |2.28| ) with the boundary conditions ( |2.30| ), it is appropriate to 
introduce the Fourier transform with respect to X\ — x 2 defined by 



Gss'in, r 2 ) 



00 H/ ~ 

—G ssl (yi,y2,iy {xi - X2) 

-00 Ztt 



(2.31) 



We then obtain ordinary differential equations in y\. Eq. (|2.29|) becomes 

d 2 



dyl 



G ss (yi,y 2 ,l) = m5(y l - y 2 )l 



(2.32) 



where k = (m 2 + I 2 ) 1 / 2 , and we look for a solution of fl2.32|) , symmetric in y± and y 2l in 
which a "reflected" part is added to the free space solution, i.e. of the form 



Zk 



(2.33) 



where A aa (I) are functions to be determined. The matrix elements G a Z „ are related to 
( [2.331 ) by the Fourier transforms of the relations ( 2.28a|) , and the Fourier transforms of 
the boundary conditions ( j2.30p determine the coefficients A aa ' (I) in ( |2.33| ). The result is 



rn 



Gll(yi,y2,l) = Yk 



GH( yi ,y 2 ,l) = - 



-iAv\-y-i\ _|_ _ 1 \ e -K(yi+y2) 
-K\yi-V2\ _ _|_ e -«(j/i+2/2) 



Gl 2 M, 1/2,/) 
Gf s {yi,y2,l) 



21 

2f 



(2.34a) 
(2.34b) 
(2.34c) 
(2.34d) 



The inverse Fourier transforms are 



GL 2 (n,r 2 ) 
Gl 2 s (n,r 2 ) 
Gf s (rx,r 2 ) 



Tfl 

— [K (mr l2 ) - K (mr* l2 ) + il(xi - x 2 ,yi + y 2 )\ 

Tfl 

— [K (mr i 2 ) - K (mr* 12 ) - il(xi - x 2 ,yi + y 2 )] 

Z7T 

in 

Z7T 

-i—I(xi - x 2 ,yx + y 2 ) 
In 



(2.35a) 
(2.35b) 
(2.35c) 
(2.35d) 



where r 12 = \t\ — r 2 |, r* 2 = \r\ — r 2 \ with r 2 = (x 2 , —y 2 ) the image of r 2 , K is a modified 
Bessel function, and / is the function (obtained as a principal value) 

- ^ sm[l( Xl - x 2 )] ^_ K(m+y2) 



J(xi - x 2 , 2/1 + 2/2) = / 

JO 



/ 



(2.36) 



Using the Fourier transforms of ( |2.28a| ) we obtain the (— s, s) matrix elements 

G\ s {y^l) = ^{- K + l + T) e ~ R{V1+V2) 
G 2 \ s { yi ,y 2) l) 



2k 1 I 



(2.37a) 
(2.37b) 



2/2,0 = ^^l-Ksign( yi -y 2 )}e-^^ + ^e-^ + y^ (2.37c) 



(</i, 1/2,0 



2k 



[-1 - k sign(yi - y 2 )] e 



-K\yi-y2\ _ 171 r -K(yi+y 2 ) 
I 



(2.37d) 



The inverse Fourier transforms are 

i(xi - x 2 ) - (yi + y 2 ) 



G_ (r\, r 2 ) 



m 
2^ 
m 



ifi(mr* 2 ) + iJ(xi - x 2 , yi + y 2 ) 



12 



<^ s , s (n, r 2 ) = — 

m 
2^ 

Gl\ s (rur 2 ) = £ 



C^ s (ri,r 2 ) 



-i(xi - x 2 ) - (yi + y 2 ) 

r* 2 

i(xi - x 2 ) - (yi - j/ 2 ) 



iTi(mri 2 ) -iJ(xi -x 2 ,y 1 + y 2 ) 



ri2 



K^mru) + iJ(xi - x 2 ,2/i + y 2 ) 



-i(xx - x 2 ) - (j/i - y 2 ) 



n 2 



K x (mr 12 ) - i J (x x -x 2 ,y 1 + y 2 ) 



(2.38a) 
(2.38b) 
(2.38c) 
(2.38d) 



where K\ is a modified Bessel function and J is the function 

J(xi - x 2 , j/i +y 2 ) = / 

Jo 



dl m MKxi-x 2 )} c _ K(m+y2) 



kI 



(2.39) 
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Since Ko(mri 2 ) diverges logarithmically as rx 2 — > 0, using ( |2.35|) in ( 2.20|) gives diver- 
gent densities, like in the free space problem. However, if we subtract from G^f its free 
space value, we get a finite result for the difference n s (y) —n s (oo) (a density depends only 
on y): 



m 



n s (y) - n s (oo) = -—K (2my) 
In 



(2.40) 



Using (|2.35|) and ( |2.38|) in (|2.23| ), we get the truncated two-body densities 



nfJ T (n,r 2 ) 



n^Zin, r 2 ) 



m 



-^-[K (mr 12 ) -K (mr* l2 )} 2 
+2J(xx - x 2 ,yi + y2)(xi - x 2 ) 



(2.41a) 



Ki(mri 2 ) KAmr 



12' 



(2.41b) 



ri2 r{ 2 

Writing (formally) the (infinite) bulk densities as n s (oo) = (m 2 /2tt)Ko(0), we see on 
( j2.4(J| ) that n s (0) = 0, as expected because of the strong repulsion between a particle 
and its image. Similarly, one checks on ( |2.41| ) that also the two-body densities vanish, as 
expected, when one of the particles is on the boundary. 

As a consequence of the assumption of perfect screening, the charge correlation func- 
tions are expected to obey a variety of sum rules. In particular, an even multipole moment 
of a particle plus its surrounding screening cloud should vanish |13[ (in the presence of 
an ideal dielectric boundary, the situation is different for odd multipole moments: some 
of their components trivially vanish for symmetry reasons, while other components have 
no a priori reason for vanishing since they are cancelled by the images). As a check of 
our results ( |2.40|) and ( |2.41| ) for the one and two-body densities, we explicitly derive the 
monopole and quadrupole sum rules. Since the one-body densities are divergent, we sub- 
tract the corresponding bulk moment for the system in infinite space, without a boundary. 
Thus, we expect the monopole sum rule 



yi>o 
n 



d 2 n 



m. ; r 2 ) - nfjZin, r 2 ) - n (2yi \r 12 ; 00) 



( 2 ) r (r* 2 ;oo) + nS(r 12 ;oo)+n^^ 



s,-sV r l2i 00 j 



-n s {y 2 ) +n s (oo) (2.42) 
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and the quadrupole sum rule 



d 2 ri y\ - (xi - x 2 ) 2 

j/i>0 
(2)T 



n 



(2)T 



(r 12 ; oo) 



... , * \ I (2)T, \ . (2)T, . 



-y 2 2 K(y 2 )-n s (oo)] (2.43) 



where n^ T (ri 2 ;oo) = — (m 2 /27r) 2 [ifo(wri2)] 2 and nf^fru', °°) = (m 2 /27r) 2 [Ki (mr i2 )] 2 
are the bulk truncated two-body densities (the moments of the system without a boundary 
have been replaced by integrals over the half-space y\ > and the addition of r^ 2 / (r* 2 ; 00) 
terms). These sum rules are derived in the Appendix. 



3 Strip 

We now consider a strip geometry. The Coulomb fluid, infinite in the x-direction, is 
constrained by two ideal dielectric walls at y = and y = W. The above studied half- 
plane case corresponds to the limit W — > 00. 

The Coulomb potential between the dielectric walls is the solution of the Poisson 
equation Av(r, r 1 ) = —27i5(r— r') with the boundary conditions for the y-component of 
the electric field E y (x,0) = and E y (x,W) = 0. By the method of images ||, after 
subtracting an irrelevant infinite constant, one gets for the Coulomb potential 

v(r, r') = — In 

where k = ir/(2W). Each particle has a self-interaction — (1/2) In | sinh k(z — z)/{ka)\. 
Writing 

sinh k(z - z') = ^ e - k{z+z,) (e 2kz - e 2kz ') (3.2) 

and using exp(2ku) and exp(2kv) instead of u and v, one can proceed as in the previous 
section, with the final result for the grand partition function, at /3 — 2, 

InS = ^Trln(l + K) (3.3) 



sinh k(z 
ka 



sinh k(z — z') 
ka 



(3.1) 
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where 



\ sinh k(z — z') sinh k(z — z') J 
In the limit W — > oo (A; — > 0), ( p~4|) reduces to ( p,18| ) as it should be. 



ka 



ka 



sinh k(z — z') sinh k(z — z') 
—ka —ka 



(3.4) 



In the continuum limit, keeping the previous definition of the rescaled fugacity m = 
27ra(/S, the eigenf unctions {^"(f)} of m _1 K and the corresponding eigenvalues 1/A are 
defined by 

- / E E ^7>y)C'V) = yCW (3-5) 
mJD z s , =±a , =l2 a 



where the domain of integration D is now the strip. Eq. 
coupled integral equations 



corresponds to a set of 



iA 



dV 



AW Jd 
-iA 



dV 



4iy 7d 
for s = ±. In terms of A, 

By using the equalities 

a. 1 



sinh k(z — z') s sinh k(z — z 

- ipi (V) + — 

sinh k{z — z') s sinh k{z — z 



lnH=-Eln 1 + X 



fl>M (3.6a) 
^ s 2 (r) (3.6b) 

(3.7) 



1 



7T 



5(r — f^) , r,r'<ED 



' sinh /c(^ — z') sinh &(;£ — z') k 
the integral equations ( |3.6|) can be transformed into the differential equations 



(3.f 



«9^ s 2 (r) 

The combination of these equations results in the Laplacian eigenvalue problem 



(3.9a) 
(3.9b) 



(-A + A 2 )V"(r) = 



with boundary conditions given by the integral equations (p76| 



(3.10) 



ti(x, 0) = -^ 2 (x, 0) and ^ s (x, W) = ^ 2 (x, W) 



(3.11) 
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For the present geometry, we look for a solution which is translationally invariant 
along the x-axis, i.e., 

^(x, y) = e ilx (A a s e Ky + B«e~ Ky ) (3.12) 

where k = (I 2 + m 2 ) 1 ^ 2 . Due to the relations ( |3.9| ), only four of eight coefficients {A", Bf} 
are independent, say {A+, 5°}. The boundary conditions (|3.11| ) imply for them a system 
of four linear homogeneous equations. The existence of a nonvanishing solution is given 
by the nullity of the determinant of this system. This gives rise to the relation between 
A and I: 



cosh 



W{l 2 + \ 2 ) 1 ' 2 



± A- 



sinh 


W{P + \ 2 ) 1 / 2 




72 + A 2)l/2 







(3.13) 



The ± sign means that, for a given /, two A's with opposite signs occur. Let us define the 
entire function 

1 



/(*) 



cosh(W7) 



cosh 



w(r + z 



2x1/2 



sinh 


w(i 2 + z 2 y/ 2 


( 


p + z 2 y /2 



(3.14) 



The solutions of ( |3.13D are the zeros of f(±z). Since /(0) = 1, we have 



m= n 

Ae/-H0) 



A 



(3.15) 



n 

ZJ-coZn Ae/ -i (0) 



m\ ( m 



Therefore, from (3.7), the grand potential per unit length to is given by 

Pu = 

= ~ [°°dl Hf(-m)f(m)\ 

Z7T JO 

In the limit W — > oo, it holds 

ln/(±m) = Vy[(/ 2 + m 2 ) 1/2 - + In [ 1 =F 



(3.16) 



in 



(p + m 2 y/ 2 / 



Infl + e-^+ofe 



-mW 



(3.17) 



Introducing a short-range repulsion cutoff l r 
the bulk term, one has finally, at f3 = 2, 



1 jo in order to avoid the divergence of 



Pu> = —f3PW + 2(3 1+ — + (e~ mW ) 
14 



(3.18) 



where 



and 



H 2vr 







fm (— ) 


+ 1 


\ma J 





m 



(3.19) 



(3.20) 



P is the bulk pressure of an infinite system ||, with a modified cutoff procedure, and 
7 is the surface tension. The leading finite-size correction term has the universal form 
7r/(24W / ). The same universal term, except for a change of sign, is found for the mass- 
less Gaussian field theory defined on the strip, with various boundary conditions of the 
conformally invariant type jn . 

We notice that the surface tension 7 can be computed directly from the density profile 
obtained in the half-space, formula ( 2.40[ ). 7 is the boundary part per unit length of the 
grand potential Q. The total number of particles is given by iV = iV + + AL = —j3(dQ/d(. 



The boundary part of this relation is 



R 9 

j3m- —7 
am 



dy [n(y) - n] 



(3.21) 



where the total particle density n(y) = n + (y) + n_(y) and n = n + (oo) + n_(oo), and we 
have used that m = 2tiQ. With respect to ( |2.40D , it holds 



[n(y) -n] = - 



m 



(3.22) 



Inserting this into (|3.21 ), one rederives the formula (|3.20|) . This result is also reproduced 
by the exact solution of the surface tension [Q (valid for an arbitrary f3 < 3 and obtained 
by using completely different means), evaluated at (5 = 2. 



4 Universality of the finite-size correction 



The finite-size correction to the grand potential, eq. (|3.18|) of the previous section, actually 
is a special case of a very general result valid at any temperature and in any dimension 
d (d > 2), for a conducting Coulomb system confined between two parallel ideal dielectric 
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plates separated by a distance W; the Coulomb system extends to infinity in the d — 1 
other directions. We shall show that the grand potential u per unit area of one plate 
(times the inverse temperature 0) has the large-P^ expansion 

Pu = —(3PW + 2/3 7 + ^ + • • • (4.1) 

where P is the bulk pressure and 7 the surface tension; these quantities are non-universal. 



However, the last term of ( Q[) is a universal finite-size correction, with a coefficient C(d) 



depending only on the dimension d: 

where T is the Gamma function and ( the Riemann zeta function. In particular, C(2) = 
tt/24. 

This universal finite-size correction is of the same nature as the ones which occur when 
the electric potential obeys periodic boundary conditions on the plates |L5] or when the 



plates are ideal conductors (Dirichlet boundary conditions) [j/j. It is remarkable that C(d) 
has the same value for ideal conductor and ideal dielectric plates. 

In the present case of ideal dielectric plates, two different derivations of (|4.1|) can 
be obtained by minor changes in the derivations which have already been made in the 
case of ideal conductor plates 0. Here, we shall concentrate on that derivation which 
relies on the assumption that the Coulomb system exhibits perfect screening properties. 
Therefore, the universal finite-size correction is not expected to hold in the absence of 
such screening properties, for instance in the low-temperature Kosterlitz-Thouless phase 
of a two-dimensional Coulomb gas. It should also be noted that if some short-range 
interactions are added to the Coulomb ones, the screening properties and therefore ( |4.1| ) 
are still expected to hold. 

For deriving ( f4.1|) from the screening properties of the Coulomb system, closely fol- 
lowing 0, we first establish a sum rule. The (i-dimensional Coulomb system is supposed 
to fill the slab < y < W, with ideal dielectric plates at y = and y = W. Let E x (0) be 
a Cartesian component parallel to the plates of the microscopic electric field produced by 
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the plasma at some point on the plate y = 0, say the origin, and let p(r) be the micro- 
scopic charge density at some point in the Coulomb system. If an external infinitesimal 
dipole p, oriented parallel to the plates, say along the x axis, is placed at the origin (on 
the Coulomb system side), since the system is assumed to have good screening properties 
it responds through the appearance of an induced charge density Sp(r), localized near the 
origin and having a dipole moment opposite to p: 

J d 2 rx5p(r) = —p (4.3) 

(choosing p parallel to the plates insures that it has to be screened by the Coulomb 
system itself, in spite of the presence of images). On the other hand, the interaction 
Hamiltonian between p and the Coulomb system is —pE x (0) and linear response theory 
gives Sp(r) = /3p(E x (0)p(r)) T ', where (. . .) T denotes a truncated two-point function of the 
unperturbed system. Thus, the correlation function obeys the sum rule 

(3 Jd 2 rx(E x (0)p(r)) T = -1 (4.4) 

We now compute the force per unit area acting on the plate y = 0. Let pd = 
(d - 2)2ix d l 2 /T(d/2) if d > 2, p 2 = 2tt. The unit of charge is defined such that the 
Coulomb interaction between two unit charges in infinite space be Vo(r— r 1 ) = \r— j J \ d ^ 2 
if d > 2 and — lnflr— f J \/a) if d — 2. Then, the Coulomb potential between the ideal 
dielectric plates can be written as a sum over images 

n=oo 

V ( y r,r')= [v (r+n2Wu- r') + v (r* + n2Wu- r')] (4.5) 

71= — CO 

where u is the unit vector of the y axis and r* = r — 2yu is an image of r. Actually, the 
sum in (|4.5|) does not converge, but it can be made finite through the subtraction of some 
(infinite) constant which we do not write explicitly since it is irrelevant in what follows. 
The force per unit area acting on the plate y = has only a component along the y axis, 
which is the yy component of the statistical average of the microscopic Maxwell stress 
tensor 

T yy (0) = -{EM - \E{rf) = - { A^(E x (0)V (4-6) 
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where we have used that, on the ideal dielectric plate, E y (0) = while all components 
of E(r) parallel to the plate give the same contribution; also, by charge symmetry, the 
average electric field vanishes and therefore (. . .) can be replaced by (. . .) T . Since the 
density n(y) vanishes on an ideal dielectric plate (because of the strong particle-image 
repulsion), the force on the plate has no contact contribution n(0)kT. Using 



4(0) 



2 dv(r, t 7 ) 



d z r 



dx' 



T> = Q 



p{r) 



(4.7) 



H) can be rewritten as 



Tyy(0) 



(d-1) 



2 dv ( r , **) 



dV 



^=0 



<4(0)p(r)> 



T 



(4.8) 



2/i d J " dx' 

As the distance W between the plates increases, dv/dx' can be expanded in powers of 
W-\ Using O in (fOJ) we find 



Tyy(fy - Tyy(fy\w=oc 



(d-l)T(d/2)((d) 



d 2 rx(E x (0)p(r)y +0 



(4.9) 



2 d TT d / 2 W d J \W d+1 . 

obeys the sum rule ([£§. Since duo/dW = T yy {0), Q gives (|Q|) 



The integral in 
and 

In the peculiar case of a one-component plasma, due to the presence of a neutralizing 
background, the bulk term in the grand potential uj is not of the form —PW. Also the 
force on a plate has an additional term related to the potential difference between the 



surface and the bulk of the plasma fT6 |, \ T7\ . However that additional term does not 
contribute to the universal finite-size correction which keeps the same form [J 



5 Concluding remarks 

The model under consideration was a two-component Coulomb gas in contact with walls 

made of ideal dielectric material. As shown in section 2, the model is solvable, in two 

dimensions at inverse temperature j3 = 2, by using the Pfaffian method. This means that, 
1 Although they have been inadvertently omitted in ref. Q, the same remarks apply to the case of ideal 
conductor plates. 
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not the grand partition function, but its square is expressible as a "treatable" determinant. 
This is the fundamental difference with the previously solved cases of the Coulomb gas 
in contact with a plain hard wall ||, or an ideal conductor wall || - |J7|, or with periodic 
boundary conditions on the plates [|1^]. As a consequence, the introduction of a four- 
component (instead of two- component) Fermi field, associated with each point of space, 
is necessary. We have worked with matrices which elements are themselves 2x2 matrices, 
so without saying it we have used a variant of the algebra of quaternion matrices. 

For the rectilinear geometry of a semi-infinite dielectric wall, we have computed the 
particle densities ( |2.40| ) and the correlation functions (|2.41| ). Due to the screening effect, 



the correlations are supposed to obey a variety of sum rules Hl3| . In the Appendix, we 
have checked the monopole ( |2.42| ) and quadrupole ( |2.43| ) sum rules, with subtraction 
of (divergent) bulk moments. The relatively complicated form of the truncated pair 
correlations prevents us, in practice, to go beyond the verification of these sum rules. 

The strip formalism in section 3 is technically very similar to the one for ideal- 
conductor boundaries 0. The main formal difference is that, when calculating the 
grand potential per unit length ( |3.16|) , an "average" (1/2) [In /(—m) +ln/(m)] instead of 
ln/(— m) should be integrated. This makes the surface tension finite, see formula ( |3.20|) . 
The universal finite-size correction term can also be obtained, as a consequence of the 
good screening properties, in the more general case of a Coulomb system of arbitrary 
dimension d (d > 2) confined in a slab of width W, at an arbitrary temperature. The 
universal finite-size correction has the same value for ideal-conductor and ideal-dielectric 
walls. In two dimensions, the ir/(24W) correction term also appears in some papers about 
the sine-Gordon theory [|18j; thus, one might suspect that the integrability of these theo- 
ries is not a necessary ingredient for obtaining this universal term, and a universal term 
might also be present in sine-Gordon models of higher dimension although they are not 
integrable. 
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Appendix 



In this Appendix we sketch the derivation of the sum rules ( |2.42|) and ( |2.43|) . Without 
loss of generality, we can choose %2 = 0. 
The l.h.s. of ( ggg) is 

— J J dyi J dxi < K (mr 12 ) K (mr* 12 ) + [Ki(mr{ 
Ki(mr '12) i^i(mr* 2 ) 



1 2 )] 2 



-2J(xi,yi + y 2 )zi 



?"12 ?T 2 

where ?/i + y 2 ), as defined by (|2.39|) , has the properties 

J{0,yi + y 2 ) = 



(A.1) 



and 



Since 



L^ e -«(?/i+y2) = m ^ ( mri2 ) 



d/ 



rax 1 



dx\ jo 
Ki(mr 12) Ki(mrl 2 



ri2 



12 



5 



(9xi 



[K (mr 12 ) - K (mr* l2 )] 



(A.2) 



(A.3) 



(A.4) 



an integration per partes on x\ transforms (|A.l ) into 

Fi = 2 | o °° d Vl jT dan {[f„K)f + [i^K)] 2 } (A.5) 

When i^o and K\ in ( |A.5| ) are replaced by their Fourier transforms with respect to x\ 

/oo I 
dl e ilxi —e~ K( - yi+y2) (A.6a) 

and 



±ixi +2/1—2/2 



Ki{mr* 12 ) 



OO 7 _|_ LT 

dl e ±itoi - Z c -«fa+?/2) 
-00 2mK 



ecomes 



27T JO J-00 

Performing first the integration on yi, we obtain 

1 



(A.6b) 



(A.7) 



m 2 roo 



m 



F 1 = — dl —e- 2 ^ = —K (2my 2 ) 

Z7T J-00 IK Z7T 



(A.8) 
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With regard to ( |2.40| ), ( |A.8| ) gives the monopole sum rule ( |2.42| ). 
The l.h.s. of ( |2~43D is 



m 2 \ 2 



)poo poo 
J dy x J dxi(yf - x\)lK Q (mr 12 ) K Q (mr* 12 ) + [K^mr 

Ki(mri2) K\{mr\ 2 



* 12 
12j 



-2J(x 1 ,y 1 + y 2 )x 1 



r 



12 



(A.9) 



Now, after an integration per partes on x±, a J-dependent term is left: 

2^ J Jo dVl J-00 dXl 1 " ^ (^°( mr i2)] 2 + [^( mr 

2 1 

H — J(x 1 ,y 1 + y 2 )x 1 [K (mr 12 ) - K (mrl 2 )] > 
m 



(A.10) 



Again, by the introduction of appropriate Fourier transforms with respect to X\ (now we 
also need the Fourier transform of xiKolmru), etc . . . ), the integral on x\ is replaced by 
an integral on /, and the integration on y\ can be performed first. After a straightforward 
but tedious calculation, the detail of which we omit, one obtains 

m 2 — I 2 m 2 — I 2 



m 2 roo 



F 2 = — I dl 
2tt Jo 



-V2 + 



-vl 



-2ny 2 



(A.11) 



An integration per partes transforms the term (m 2 — l 2 )y 2 / into 2l 2 y\j ^ . Thus 



, m 2 r°° dl 



,m 



F2 = yl^ / - e" 2 ^ = y 2 2 ^-K (2my 2 ) 

Z7C JO K 27T 



(A.12) 



With regard to (|2.4(J|), (|A.12|) gives the quadrupole sum rule (|2.43|). 
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